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Abstract. The new hadronic currents implemented in the TAUOLA library are obtained in the 
i-£h I unified and consistent framework of Resonance Chiral Theory: a Lagrangian approach in which 

the resonances exchanged in the hadronic tau decays are active degrees of freedom included in a 
way that reproduces the low-energy results of Chiral Perturbation Theory. The short-distance QCD 
, constraints on the imaginary part of the spin-one correlators yield relations among the couplings 

—< that render the theory predictive. 

In this communication, the obtaining of the two- and three-meson form factors is sketched. One 
of the criticisms to our framework is that the error may be as large as 1/3, since it is a realization 
of the large-Afc limit of QCD in a meson theory. A number of arguments are given which disfavor 
(yT) \ that claim pointing to smaller errors, which would explain the phenomenological success of our 

description in these decays. Finally, other minor sources of error and current improvements of the 
" code are discussed. 
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INTRODUCTION 

New hadronic form factors have been included in TAUOLA [1], the standard Monte 
Carlo generator for tau lepton decays [2]. In this contribution we concentrate on the 
theoretical inputs and the associated errors. Practical aspects that may be interesting for 
the user are the topic of O. Shekhovtsova's communication [3]. This project is essential 
to meet the experimental requirements, as discussed in Ref. [4]. The definition we use 
for the hadronic form factors, J^, in the decay to a hadronic state H, x~ (P) — > Hv r (N) , 
is 

V2 

We will consider the two- and three-meson hadronic states in turn 



Speaker. 



p. - -^(A0/(l-r 5 )"(P)^. CD 
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FORM FACTORS IN TWO-MESON x DECAYS 



For T — > [hi {p\)h2{pij\ v? decays, the hadronic current reads 



(pi -VI- ^r(P! + V(*) + ^( pi+P2 )»F p (s) 



(2) 



where s = (pi + P2) 2 and A12 = m\ — m\, with N K K ° = 1. The other final states are 
given by ST/ (3) symmetry. The dominant contribution to these decays is given by 
the vector form factor, F v (s), which is obtained using Resonance Chiral Theory [5] 
(RXT). It reproduces the low-energy results of Chiral Perturbation Theory [6] (%PT) 
at NLO in the chiral expansion and includes the light-flavored resonances as active 
degrees of freedom in the theory without any ad-hoc dynamical assumption. R%T is 
a realization of the large-Afc limit of QCD [7] in a theory for the lightest mesons and 
resonances. When the asymptotic vanishing of the form factor is imposed on the result, 
F v (s) = \+F v G v /F 2 -s/(M^-s), it yields F v (s) = M^/(M 2 -s) = F VMD (s), the 
vector meson dominance result. 

FSI among the two pseudoscalar mesons are encoded in the %PT loop function, 
Apq(s), whose imaginary part enters the vector meson off-shell width [8] included via 
F VMD (s) = My /(My — s — iMyTy(s)). There are several ways to resum them, such 
as the exponentiation of the real part of the Omnes [9] function [10] or the use of 
dispersion relations [11, 12]. The approach proposed in Ref. [13] and used in Ref. [1], 
takes Ref. [10] as a starting point and includes phenomenologically the effect of the 
excited p-like resonances when they are seen in the data, as in the two-pion mode [14]. 
The general structure for the single-resonance contribution is 



F v P0 {s)=F VMD {s)^ V _ 



(3) 



where Y*p,q extends over the light pseudoscalars with suitable quantum numbers to 

contribute in the loop and N^® is dictated by SU (3) symmetry, with Nf 00 ^ = 1. 
The scalar form factor is important to describe correctly the low-energy region of the 
data [15] for the T~ — > (Kn)~v T decay [16]. In the 2012 release of our code distribution 
for [1] it has been included following the coupled channel analysis of Refs. [17] for 
the strangeness changing scalar form factors. The strangeness conserving scalar form 
factors [18] may be important in the modes including an r\ -meson. 

FORM FACTORS IN THREE-MESON x DECAYS 

In the T~ — v [hi{pi)h2(p2)hT,(p'i)\^V x decays, the hadronic current reads 

+ c ^ Mi7 4 - ^^2C 5 e% p v lP p 2 p^F 5 1 , (4) 
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where q^ 1 = (p\ + P2 + P3)* 1 , F is the pion decay constant and N = (Vfj ) /F. Only 
three among the F\ , F 2 and F 3 form factors are independent. The choice for each mode, 
with the Ci factors, can be found in Ref. [1]. 

In this reference, the 3k and KKiz modes were included within the framework of R%T 
following Refs. [19] (where T ai (q 2 ) can be found) and [20], respectively (the mode 
r\K~K Q [21] will be available soon). In these decays, vertices with more than one 
resonance are considered, including the VVP vertices [22] (P stands for one of the 
lightest pseudoscalar mesons) and the VAP vertices [23]. Some of the abundant new 
couplings can be related upon imposing that the imaginary part of the VV and AA 
correlators go to zero asymptotically in the large-Afc limit [24]. It is noteworthy that 
the relations found are in agreement for all three cases, and with those obtained in the 
one-meson radiative tau decays [25]. 



DISCUSSION ON THE ERRORS: AS LARGE AS 1/3? 

The smallness of a expansion parameter for applying perturbation theory is given by 
the size of the coefficients of the expansion. Usually, in perturbation theory, i.e. QED 
or %PT, it is easy to compare the LO and NLO expressions to determine the expansion 
parameter and the smallness of the related coefficients warranties a good convergence 
of the lowest order computations to the true result. It is not possible to resum the infinite 
number of diagrams that appear in the Nq — > 00 limit of QCD to be able to judge if the 
expansion parameter is small enough to rely on this approach. One can only derive that 
at NLO the non-planar diagrams are suppressed as \/Nq and the diagrams with internal 
quark loops as l/Nc- Moreover, a factor of nf could enhance the latter diagrams, but 
the fact that there is negligible mixing between the qq and the qqqq states hints that 
these kind of diagrams are heavily suppressed by their coefficients. Therefore, I/Nq 
would be a better estimate of the effective expansion parameter, which agrees with the 
phenomenological success of its predictions on the hadron side. In particular, we note 
the good convergence of the predictions of the %PT couplings working at NLO in the 
l/Nc expansion within R%T [26] and the successful description of the hadronic decays 
of the tau lepton as indications that the expansion parameter is indeed smaller than 1/3. 
Noticeably, the actual expansion parameter can be computed for R%T in the study of the 
vector form factor of the pion at NLO in the l/Nc expansion [27], yielding 

2 F 2 96nF 2 ' K } 

which, at lowest order, is the ratio of the vector width and mass, ~ 0.2, agreeing with 
the previous discussion. Moreover, we should emphasize that our approach goes beyond 
the Nc — >• 00 limit. We supplement the lowest order in the l/Nc expansion for the theory 
in terms of mesons by the leading higher-order correction, namely by including the 
resonance (off-shell) widths for the wide states p, K* and a\. This seems to point to 
smaller errors (< 10%) than those characteristic of the LO contribution in the l/Nc 
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expansion and may be able to explain, altogether, fine agreement with the data 2 . 



Other sources of error 

The procedure outlined in eq.(3) [10], violates the unitarity and analiticity constraints 
at NNLO in %PT . However, it is possible to devise a strategy which overcomes this 
problem, as put forward in Ref. [12]. Within it, both the real and imaginary parts of the 
%PT loop functions are kept in the denominator of the form factor and the tangent of 
the relevant phaseshift, 8j(s), is defined as the ratio of its imaginary and real parts and 
it enters a three- subtracted dispersion relation in order to provide the final form factor. 
For the K% [12, 16, 28] and nn modes [1 1, 13, 29] the numerical differences are smaller 
than the experimental errors. For the decay modes where the elastic approximation is not 
good (such as Kr/W), eq.(3) can be an approximation which induces a small error while 
a coupled channel analysis is developed. For the Kk, nn and ltr\ W - which proceeds 
dominantly through the it — r\ — TjW mixing [30]- both approaches could be employed 
and the error induced by using eq.(3), coded in the TAUOLA 201 1 version [1], is at the 
percent level. 

BaBar data for the T~ — > n~n~n + v T decay [31] has confirmed the nice description of 
Ref. [19] for dT/dq 2 improving the performance of old TAUOLA hadronic currents (it 
happens similarly in the two-meson modes) [32] 3 . However, a deviation in the low- 
energy region of the differential decay width as a function of the Tl + Tl~ invariant mass 
remains. It is possibly due to the absence of the a meson in the theoretical description. 
Its thorough incorporation requires the inclusion of inhomogeneities as angular averages 
of the form factors [34] but this is very time consuming and not practical for the Monte 
Carlo. Instead, we include [32] an educated parametrization of the energy dependence 
of the / = 0, 2 contributions [35] and use 8q(s) and 62(5) consistent with the chiral con- 
straints at low energies [36]. The reduction of the errors in that region will be presented 
in Ref. [32]. 

Other sources of error are smaller, such as: the contribution of SU(3) breaking terms 
[37]; the inclusion of excited resonances from a Lagrangian [38, 39], and, in the %% 
mode -where data are more precise- SU(2) breaking [40]; and some NNLO [41] sub- 
dominant terms which are not reproduced [42] by Guerrero-Pich-like parametrizations 
[10]. Their effects are discussed in more detail in Ref. [1]. Some improvements, along 
the lines discussed above will be included in the next release of our currents, see 
Ref. [32] for details. 
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It remains exciting to confront this well founded conjecture with experimental data in new applications. 
3 MC-Tester [33] is a useful tool for this kind of comparisons. 
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